Phase-space behavior and conditional dynamics of an optomechanical system 
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We characterize the local properties of an optomechanical system comprising the movable mirror 
of a resonator and its intracavity field, mutually coupled via radiation-pressure. Our approach 
shows that both the state of the mirror and the field can be interpreted as squeezed thermal states 
whose dynamical properties can be tuned by properly choosing the working parameters. This 
allows us to design conditional procedures for the amplification of the correlation properties of the 
optomechanical state. Our study is a step forward in the understanding of the physics that rules a 
system of current enormous experimental interest. 
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The observation of quantum dynamics of a macro- 
scopic system is an outstanding goal in modern physics. 
Several paths have been followed in order to achieve 
macroscopic quantum mechanical behaviors, using nano- 
electro-mechanical devices coupled to single-electron 
transistors or optomechanical systems coupling electro- 
magnetic fields and vibrating mechanical structures [l]. 
The desideratum common to all such physical situations 
is the achievement of negligible thermal excitations of any 
mechanical system, whose energy has to be pushed down 
to a level comparable with their specific single-quantum 
transition energy. Recently, some notable results have 
been achieved in this respect, both theoretically and 
experimentally [2|-[6|, by exploiting back-action induced 
mechanisms. Moreover, correlations of quantum mechan- 
ical nature have been predicted in micro-optomechanical 
systems, resilient to the exposure to Markovian envi- 
ronmental noise and relatively high operating temper- 
atures [2]. Progresses towards the experimental revela- 
tion of their classical counterparts have been recently re- 
ported [7] . This complements previously studied unitary 
dynamics [9] and sets the ground for realistic experimen- 
tal investigations. Indeed, operative schemes for the ob- 
servation of optomechanical entanglement have been sug- 
gested, based on the use of ancillary fields and all-optical 
observations [10]. Noticeably, very recently, evidences of 
quantum effects in the motion of a mechanical system 
have been observed [8[. 

In the optomechanical setting comprising an optical 
cavity with a movable mirror coupled to light through ra- 
diation pressure forces, however, a complete understand- 
ing of the effects that the coupling with light has on a 
micromechanical device is still lacking. The "not im- 
mediately intuitive" nature of the correlations studied 
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in [H, [lo|, for instance, deserves a deeper investigation: 
the properties of the whole system and its constituents 
have to be properly unravelled. In this context, in con- 
sideration of the formal equivalence between the movable 
mirror and a quantum harmonic oscillator, it is natural 
to use techniques and tools of quantum optics to provide 
a more explicit interpretation of the dynamics of the sys- 
tem. This is exactly the direction considered in this pa- 
per, where we draw a phase-space analysis of the cavity 
mirror-field system depicted above. This allows not only 
for a complete characterization of both the mirror and 
the field as squeezed thermal states, but also the identifi- 
cation of conditioning strategies, based on measurement 
and postselection, through which the state of the mirror 
is significantly modified. 

The paper is organized as follows. In Sec. |T] we intro- 
duce the formalism and technical details being used in 
order to tackle the problem at hand. A detailed and 
extensive description of the methods can be found in 
Refs. [i,[io| and here we only sketch the strategy to follow 
in order to gather the form of the covariance matrix of 
each optomechanical system. Seclllis devoted to the de- 
scription of the dynamics, in phase space, of the reduced 
state of the mirror and the field of one of the cavities. We 
determine occupation number and degree of squeezing of 
the equivalent squeezed thermal states used in order to 
effectively describe the reduced state of each subsystem. 
In Sec. Hm we track the changes undergone by these fig- 
ures of merit when the optical part of the system or a 
proper ancilla is subject to conditioning measurements 
such as homodyne detection or projection onto the vac- 
uum. Finally, Sec. HVl summarizes our results. 



I. INTRODUCTION TO THE FORMALISM 

We consider an optical cavity with a movable, highly 
reflective end- mirror, coupled to the intracavity field by 
radiation pressure [11]. The Hamiltonian of the system 
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reads (we take units such that h = 1 across the paper) 
n = {uJc-Goq)a^a^'^{p^^f)-i£{ae'''^'-h.c.), (1) 

where a and are the cavity-field bosonic operators, Uc 
its frequency and ujm the frequency of the mirror, which 
is modelled as a harmonic oscillator with quadratures 
q = {m^ + rh)/^/2 and p = i{rh) — rh)/\/2 and bosonic 
operators m, rh) . Moreover, 8 is the coupling rate be- 
tween an external driving field (with frequency uOq — ^c) 
and the cavity field and Go is the radiation-pressure cou- 
pling strength. Starting from ?i, Refs. 0, |5|, [lOj show 
how to obtain a set of linear Langevin equations useful 
to the reconstruction of the dynamics of the optomechan- 
ical system. This can be done by focusing our attention 
on the fluctuations of the relevant operators in the prob- 
lem. The dynamics of the system is formally described 
by the integral matrix equation 
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where f(t) = {5q Sp Sx Sy)^ is the vector of the quadra- 
ture fluctuations and n{t) = (0 ^{t) V2kSx''' V^ilSy'"')^ 
accounts for the noise entering the system. Here, ^(t) 
describes the Brownian motion of the mirror due to its 
coupling with a phononic bath at temperature T and 
fe^^ = ((5a^^t + 3a''') /V2 and 5y''' = i(^a^^t _ Sa''')/V2 
are the quadrature fluctuations of the input noise to the 
cavity, which has an amplitude decay rate k,. In the or- 
dered basis given by {Sq^ Sp^ 5x^ Sy}^ matrix K reads 
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where G = G^EI\J + is a modified mirror-field cou- 
pling rate and A is the radiation-pressure affected sys- 
tem detuning [H, [5| . In the parameter regime where the 
eigenvalues of K have negative real parts, Eq. (j2j) can 
be solved at the steady state defined by lim^^oo = 
0. The linear nature of the problem at hand guar- 
antees the Gaussian-preserving character of the map 
defined by Eq. ([2]). This allows us to focus on the 
entries of the covariance matrix V defined as Vij = 
limt^oo ^{{fi{t)^ fj{^)}) [l^l, where the expectation value 
is calculated over the mirror-field state. The explicit form 
of V is 



V = 
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where M [F) is the 2x2 block matrix describing the lo- 
cal properties of the mirror (field) and C accounts for the 
mirror-field correlations. By assuming a large mechanical 
quality factor and a Ohmic spectral density of the back- 
ground phononic bath [2i, |lo|, Hll , analytic expressions of 



the elements entering V can be easily found. They are 
however too lengthy to be reported here. As discussed 
above, Eq. (|2]) does not affect the Gaussian nature of any 
input mirror-field state. In Refs. [5,^0], this property has 
been used to show the entangling properties of the radia- 
tion pressure coupling in the presence of noise. Here, we 
concentrate on the phase-space properties of the reduced 
state of each subsystem. 



II. RECONSTRUCTION OF THE WIGNER 
FUNCTIONS 

In this Section we give an effective description of the 
local states of the mirror and the field in terms of equiva- 
lent linear optics transformations applied to proper input 
Gaussian states. An analogous study performed with re- 
spect to the joint two-mode state described by Eq. (H]) 
is made difficult by the fact that V describes a mixed 
state. Finding the corresponding linear-optics interfer- 
ometer from which the state would be derived is, in gen- 
eral, a demanding task which goes beyond the scope of 
our investigation. A complete and self-contained account 
of the quantum correlations within V has already been 
provided [2]. Here we concentrate on the reduced state 
of the mirror and the field as resulting from the radiation 
pressure-induced evolution. The inspection of Eq. (j4]) is 
quite revealing. Indeed, by looking at the form of M, 
one notices that the mirror could be in a squeezed state 
with squeezing factor | ln(^^). In the phase space, such 
squeezing would necessarily fee along the direction of one 
of the two mirror quadratures {i.e. the squeezing param- 
eter must be real). On the other hand, the state of the 
field could be squeezed as well, but this has to occur along 
a direction that can be completely general. A formal ac- 
count of such effects is given in Appendix A. Here, we 
make our expectations quantitative and concentrate on 
the behavior of the Wigner function associated with each 
reduced subsystem. These are related to the covariance 
matrices introduced in Eq. (j4]) by the expression 
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where J = M or F. Here, /m = {^q ^p), fr = {^x Sy) are 
complex vectors for the mirror and field quadrature vari- 
ables. For the mirror, we have taken Sq = {a -\- a*)/>/2 
and Sp = i{a* — a)/V^ with a = ar -\- ica and the field 
quadrature variables Sx and Sy have been analogously 
defined. With these choices, WM(<^r,<^i) turns out to be 
a function of the effective detuning A. In Fig.[T]we show 
such a dependence by setting A and studying WM(<^r, <^i) 
against ar and ai. For a resonant interaction [panel (a)], 
WM{(^r^(^i) is an isotropic Gaussian with a large width 
(notice the range of values of a^^^), which is indicative 
of a thermal state with large mean occupation number. 
As A grows, the variance of WM(<^r, <^i) along the ar di- 
rection becomes sensibly larger than that along a^, thus 
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FIG. 1: Wigner function associated with the mirror state plot- 
ted against the rescaled detuning A = jujm/lO {j G Z). Each 
panel corresponds to a set value of j. We have taken j = (a), 
4 (b), 7 (c), 10 (d), 13 (e), 16 (f), 19 (g) and 20 (h). Through- 
out the paper we have used the following set of parameters: 
cj^/27r = lOMHz, ujo - uJc/27t = 3.7xlO^^Hz, 8 = GxlO^^^z^ 
Go 1400Hz, T = 0.4K, and a cavity of finesse 10^ that is 
1mm long. 



showing a A— dependent squeezing effect. However, the 
direction of squeezing remains fixed, thus indicating that 
the squeezing parameter is always real and negative. The 
width of the Wigner function, which gives information on 
the effective temperature of the mirror, is considerably 
reduced in a large range of values of A, thus revealing a 
strong cooling effect, even for moderate detuning between 
the cavity and the external field. This effect has been 
extensively studied in a series of recent experiments @. 
For large values of A, the state of the mirror goes back 
to a thermal behavior, with a reduced mean occupation 
number, which slowly increases toward its initial value 
as A is further increased. In fact, with a growing de- 
tuning, the effective coupling rate between mirror and 
field decreases as less input power couples into the far 
off-resonant cavity, thus driving the mirror towards its 
unperturbed initial thermal state. 

In order to make these considerations more quantita- 
tive, we proceed as follows. The most general single- mode 
Gaussian state of a boson of frequency uj is given by the 
squeezed thermal state 

Pgen = ^(/?)^(re^^)e-|(2n+i)5t(^e«v), (6) 



where Z{P) = e^l^ — e is the partition function of a 
thermal state having effective temperature = ki)T /uj 
{kb is the Boltzmann constant), S{re^'^) is the single- 
mode squeezing operator with complex amplitude re*"^ 
defined in Appendix A and ft is the bosonic number op- 
erator. Here, we have neglected the possibility of a dis- 
placement of the state in the phase-space as the quadra- 
ture operators we are considering refer to the zero-mean 
fluctuations of both the field and mirror. In order to 
quantify the similarity between the mirror state and pgen , 
we could compare the behavior of the associated Wigner 
functions. However, as we are dealing with Gaussian 
states, a much more convenient comparison is performed 
in terms of covariance matrices [l^. We use sje*'^^ and 
nj + 1/2 to indicate the effective squeezing factor and 
thermal variance of mode J = M, F. In this way, a full 
characterization of each reduced state is obtained by nu- 
merically solving the matrix equation 

J = vi = S{sj, ifj){nj + ^)l25(sj, (7) 

for the unknown effective parameters nj sj and 
where S{sj^ipj) = {s^ cp) is the symplectic transforma- 
tion corresponding to single-mode squeezing of amplitude 
sje^"^-^ [cfr. Eq. ()A-ip ] and (nj + ^)l2 is the covariance 
matrix of a single- mode thermal state (I2 is the 2x2 
unit matrix). Any explicit dependence of nj and sj on 
A has been omitted for convenience of notation. The so- 
lution of Eq. ([7]) for J = M leads to the behaviors shown 
in Fig. O As expected, we find ipM = tt, in line with our 
conjecture on the antisqueezing of the mirror state. On 
the other hand, um and sm are not monotonous against 
A. In particular, the effective mean occupation number 
is peaked in correspondence with the maximum of sm- 
Although our numerical approach allows for the immedi- 
ate visualization of the results, it is possible to give an 
analytic expression for both the degree of squeezing and 
the mean occupation number. We refer to Appendix B 
for full analytical details. 




FIG. 2: Effective squeezing factor sm (rescaled so as to im- 
prove its visibility) and average phonon number tTm against 
the dimensionless detuning A/ujm. The inset shows tTm with 
an enlarged vertical axis so as to provide evidence of the con- 
siderable cooling experienced by the mirror. We find cpM = vr, 
regardless of A. 
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A natural question arises, at this point: How is the 
effective degree of squeezing related to the optomechan- 
ical entanglement that is known to occur in the condi- 
tions at hand? In order to provide a quantitative an- 
swer, we have calculated the logarithmic negativity of 
the mirror- field system [i,[io| as a function of e for 
sm going from to its maximum 1.1 and back to zero 
again. Two different values of entanglement will be asso- 
ciated with the same degree of squeezing 5m, depending 
on the way the detuning has been changed. This is an 
effect of the asymmetry, with respect to its maximum, of 
the curve showing entanglement against detuning [Fig. [3] 
(a)]. Such considerations are confirmed by the incom- 
plete hysteresis loop-like curve displayed in Fig. [3] (b). 
In each branch of the hysteresis loop, entanglement is 
a monotonic (increasing or decreasing) function of the 
squeezing. The qualitative interpretation of our results 
goes along the following lines. In the linearized regime 
where Eq. ([2|) is valid, the radiation-pressure coupling 
term —God^aq in H includes both energy-preserving and 
non-preserving terms [15]. On their own, the energy pre- 
serving terms, having the form difi^ h.c.^ would gener- 
ate a beam-splitting operation involving mirror and cav- 
ity field. The non-preserving term b)rh) + h.c. is the 
generator of a two- mode squeezing transformation (23| . 
Under this viewpoint, the full optomechanical dynam- 
ics encompassed by Eq. ([1]) and without the inclusion of 
noise effects, can be interpreted in terms of an equiva- 
lent interferometric setting where two modes impinge at 
a beam splitter. It is well-known that, in order for the 
output of such a device to show entanglement, a certain 
degree of non-classicality (such as squeezing) should be 
present at the input. Therefore, maximum entanglement 
is achieved at the conditions that, dynamically, optimize 
the squeezing of the mirror mode. Clearly, the true phys- 
ical situation at hand is somehow different from this in- 
tuitive picture as we deal with a stationary state that 
is affected by losses and noise. The above discussion, 
however, is sufficient to get an intuition of the physical 
process at hand. 

Let us now shift our attention to the cavity field state 



(a) (b) 




FIG. 3: (a): Modulus of the effective squeezing factor sm 
and entanglement in the mirror- cavity system plotted against 
/\/ujm- The two curves are peaked at the same detuning, 
(b): Logarithmic negativity in the mirror-field system 0, [13] 
shown against e~^*^. The horizontal axis is truncated at 
Sm — 0.3 for clarity. 



and consider the Wigner function ^^^(a^,^^) associated 
with its reduced state. This is easily calculated by means 
of Eq. with J = F and its behavior is shown in Figs.H] 
for set values of A. It is interesting to note that squeezing 
appears in the field state even at zero detuning, as wit- 
nessed by the evident deformation of the bidimensional 
Gaussian in Fig. |4] (a). This is in contrast with the be- 
havior of the mirror state, which requires A 7^ in order 
to be squeezed. Therefore, as it will be clarified later, 
one cannot relate the squeezing of the field to the entan- 
glement in the joint cavity-mirror state. The cavity field 
squeezing appears to be a simple consequence of the cou- 
pling to the mirror. Indeed, we have studied the evolution 
of the field squeezing at A = and by replacing the cou- 
pling parameter G with g = xG, where x ^ i^A]- This 
accounts for a resonant cavity-pump configuration and a 
tunable power of the cavity-driving field. The shape of 
WF{(^r^(^i) for proper values of x is shown in Figs. O 
where we see the rapid appearance of squeezing (in panel 
(c) we had to enlarge the range of ai up to [—500, 500] 
in order for the Wigner function to look isotropic), more 
pronounced as the scaling factor increases. The sudden 
squeezing of the cavity field state, which is absent in the 
mirror state, may be due to the different initial condi- 
tions assumed for the two subsystems. The explicitly 
mixed nature of the mirror state (due to the thermal 
background of phononic modes) requires the enhance- 
ment of the nonlinear character of the interaction with 
the field in order to exhibit significant squeezing. This 
is effectively achieved via a non-zero detuning. On the 
other hand, as the cavity field is prepared in a coherent 
state, even a resonant interaction significantly affects its 
state. 

By following the same line as in the case of the mirror, 
we can calculate the effective mean occupation number 
and degree of squeezing of the field state against the op- 
tomechanical coupling and the detuning, so as to gather 
a complete picture (see Fig. [6|). It is easy to show that 



uf = h V det F, sf = -arccosh f — , ) , 

2 2 \2VdetFj 

(8) 

which can be easily calculated by using Eq. dH and the 
approaches of Refs. P, IToj. Any resulting expression is 
however too lengthy to be reported here and we summa- 
rize our findings in Fig. [6l First, ni? at A = appears to 
be a quadratically-increasing function of x^ as we have 
checked by finding the best fit function corresponding to 
this configuration [cfr. panel (a)]. The effective squeez- 
ing, on the other hand, starts from zero as a concave func- 
tion, soon becomes convex and grows as y^, as shown in 
Fig. [6] (b). This proves that even at A = squeezing of 
the field state should be expected. On the other hand, as 
previously commented, we have checked that no squeez- 
ing appears in the mirror state at resonance. At A = 0, 
tuning X simply determines the isotropic increase of the 

width of WM{(^r^ (^i)' 

We now address a second interesting point: different 
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FIG. 4: Wigner function associated with the cav- 
ity field plotted against A. He have taken j — 
0, 0.4, 0.7, 1, 1.3, 1.6, 1.9, 2 in going from (a) to (h), 
respectively. 



FIG. 5: Wigner function associated with the cav- 
ity field studied against the parameter x that tunes 
the optomechanical coupling. We have taken x — 
0.01, 0.04, 0.4, 0.5, 0.6, 0.8, 0.9, 1 in going from (a) to (h), 
respectively. 



from the case of the mirror state, at x = 1 

rotates (in the phase space) with the degree of squeezing, 
which implies that the squeezing factor of the field state 
is complex and provides additional information on the 
mechanism for setting entanglement within the system. 
In fact, by using again the equivalent interferometric de- 
scription depicted before, one can conclude that along 
with the amount of squeezing of the mirror state, an im- 
portant role is played by the relative direction of squeez- 
ing of the two subsystems states. While the state of the 
mirror is squeezed along a fixed direction in phase space, 
the direction of squeezing of the field changes with the 
detuning as shown in Fig. SI thus affecting the amount 
of entanglement in the output modes. In this picture, 
we expect the detuning to enter in a critical way in the 
splitting ratio of the effective beam splitter that superim- 
poses the two modes. It is important to note that a study 
including larger values of x is not possible unless the set 
of parameters we are using is adjusted. Indeed, larger 
values of x soon put the system in a regime of instability 
where a stationary solution to the Langevin equations is 
no longer possible. As we would like to keep the level 
of technicalities related to the search for stability regions 
away from our discussion, we refer to Refs. [2|, [lOj] for a 
detailed description of this issue. 



III. CONDITIONAL MEASUREMENTS 

We now address the second relevant point of our study, 
namely the search for measurements on the accessible 
part of the system that allow to condition the state of 
the mirror. For instance, by taking the degree of squeez- 
ing of the reduced state of the mirror as an indicator 
of its non-classicality, is there a way to push it closer 
to the vacuum limit? Here we consider measurements 
performed over the field interacting with the mechani- 
cal mirror as well as an ancillary mode that is superim- 
posed to the field state at a beam splitter. One can in 




FIG. 6: Effective variance (a) and degree of squeezing (b) of 
the equivalent field state plotted against A/cJrn and x- 
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fact find the covariance matrix of tiie system compris- 
ing tiie mirror and ttie extra-cavity field. It is related to 
the intra-cavity field studied so far via the well-known 
Collett-Gardiner input-output relations [16]. A detailed 
calculation (see Ref. [10]) reveals that, with a proper 
temporal/frequency- filtering of the quadrature values of 
the optomechanical system, V in Eq. (|4]) can be recon- 
structed from the (experimentally determined) output 
covariance matrix Vout as F ex (Fo^^ — 1/2). In Ref. [17] it 
was shown that this technique can also be used to have ro- 
bust entanglement between the output mode and the me- 
chanical mirror. More recently, the pulsed optomechan- 
ical paradigm has been proposed as an effective mean 
to reconstruct the state of a mechanical device [l8|. A 
first proof-of-principle demonstration has been reported 
recently [19]. 

In any case, the linearity of the relations at hand al- 
lows to work directly with V. We restrict our attention 
to Gaussian measurements, which preserve the Gaussian 
character of the state being measured. We explicitly con- 
sider projective measurements of a field mode onto the 
vacuum and homodyne measurements. The changes at 
the level of second moments of the quadrature operators 
after either a projection onto the vacuum or a homo- 
dyne measurement can be given in terms of Schur com- 
plements [20]. Let us assume we have the following co- 
variance matrix of ni + n2 modes 

Ve. = ^) , (9) 

where B gives the local variances of quadrature operators 
in a system of n2 modes to be measured, A describes 
the analogous quantities for the remaining ni modes and 
C accounts for the correlations between the ni- and n2- 
mode systems. For a homodyne measurement over the n2 
modes, the "updated" covariance matrix of the remaining 
rii modes {i.e. the covariance matrix of the system after 
the measurement has been performed) is given by the 
Schur complement 

AL™ = A - C[7rB7r]-i^,C^, (10) 

where 7r = ^'J^-^ ^ ^ and [•]~p^ stands for the Moore- 

Penrose pseudo-inverse of a matrix [l2|, [20| . On the other 
hand, if the 712 modes are projected onto the vacuum, the 
updated covariance matrix of the remaining subsystem is 

a;„, = a-c(b + i)-ic^, (11) 

where 1 is the 2n2 x 2n2 identity matrix and the stan- 
dard inverse of a matrix is used. While C is a 2n2 x 2n2 
matrix, both ttBtt and B + 1 have dimension 2ni x 2ni. 
The term to be subtracted to the pre-measurement co- 
variance matrix A is always a 2ni x 2ni matrix, as it 
should be. We first consider a measurement to be oper- 
ated directly on the cavity field. For the case at hand, 
ni = 77-2 = 1, A = M, B = F and the calculation of the 



updated covariance matrix of the mirror is quite straight- 
forward [21]. Through an equation analogous to Eq. ([7]) 
but involving, this time, M^^^ we easily get infor- 
mation about the effective mean occupation number and 
squeezing of the conditional state of the mirror. 

The comparison between the case of vacuum-projected 
(solid line), homodyne- measured (dashed line) and non- 
conditioned (dot-dashed line) mirror state is shown in 
Fig. [7] (a). Effective cooling of the mirror is achieved, 
for both the homodyne measurement and the projec- 
tion onto vacuum. At A O.ScJm, a reduction of um 
of almost one order of magnitude with respect to the un- 
conditioned case is achieved when a vacuum-projection is 
implemented. This is accompanied by the raising of the 
effective squeezing parameter at small detuning, as shown 
in Fig. [71 (b). In this sense the projection onto vacuum, 
being far from a "classical" measurement such as homo- 
dyning, pushes the mirror state towards non-classicality. 
Cooling and squeezing effects are both accompanied by a 
rotation of the Wigner function associated with the state 
of the mirror (see Fig. [8l where the large initial squeezing 
and its decrease to zero are clearly evident). 

Although effective in pushing the mirror towards a 
more non-classical state, the strategy used above obvi- 
ously cannot be used to affect the optomechanical en- 
tanglement. In order to achieve this task, we have to 
consider an ancillary mode that effectively introduces an 
additional control. We thus consider a second field, su- 
perimposed to the one studied so far at a beam splitter 
with a tunable splitting ratio 0. The ancilla is initially 
prepared in its vacuum state and, consistently with the 
analysis above, projections onto vacuum are considered. 
The mixture of the cavity field with the ancilla can be 
described formally by considering the analogy between 
injecting a field into a cavity and a beam splitting oper- 
ation Bs = ex.-p[6{a^b — /i.e.)], where b is the annihilation 
operator of the ancilla and is the splitting ratio. 

In order to correctly understand what has to be ex- 
pected from this thought-experiment, we have to analyze 
the dynamics of the entanglement within the tripartite 
system comprising the mirror and the two fields. This is 
easily done by considering the covariance matrix of the 



(a) (b) 




FIG. 7: Panel (a) shows the effective mean occupation num- 
ber of the mirror state against A/cJm- Panel (b) reports 
the effective squeezing. Dot-dashed lines show the proper- 
ties of the unconditioned mirror state, dashed lines are for 
the homodyne-measurement case and solid lines are for the 
vacuum-projection case. 
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FIG. 8: Snapshots of the evolution of the Wigner function as- 
sociated with the mirror state against the detuning A [which 
increases from to lOcJm in going from (a) to (h)] for the 
case of projection onto vacuum of the cavity field. 



system after the beam splitting operation superimposing 
the cavity field and the ancilla and evaluating the loga- 
rithmic negativity of each of bipartition extracted from 
the tripartite system. The analysis can be performed 
analytically to some extent and, in principle, complete 
expressions for the logarithmic negativity of every re- 
duced bipartite state can be given, although they are 
very complicated. We find that no entanglement can be 
established between the cavity field and the ancilla, re- 
gardless of the splitting ratio used for the beam split- 
ter and the detuning in the cavity-mirror system. The 
two field modes remain separable. On the other hand, 
there is an interplay between the entanglement in the 
mirror-cavity field and in the mirror-ancilla subsystem. 
In general, both can be entangled, thus revealing a one- 
mode biseparable (or two-way entangled) nature of the 
trimodal state we are analysing, where the classification 
given in [22^ is used. This is somehow complementary 
to the three-body scheme that has been considered in 
Ref. [io|, where a genuine tripartite entangled state is 
created. The quantitative results relative to the two-way 
entangled state we get here are shown^in Fig. O where 
we present the logarithmic negativity [2^ between mirror 
and cavity field (dashed line) and between mirror and an- 
cillary field (solid line) against the detuning, for different 
values of the splitting ratio 0. The entanglement ini- 
tially present in the cavity-mirror system is poured into 
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FIG. 9: Entanglement in the mirror-cavity (dashed line) and 
mirror-ancilla system (solid line) against the detuning param- 
eter and as is changed from to ^ tt. Entanglement is cre- 
ated in one subsystem at the expenses of the other one. In 
panel c ((e)), ^ :^ 7r/4 (^ ^ 37r/4), while panel (d) shows the 
case of close to 7r/2. 



the ancilla- mirror one as goes from to 7r/2, when 
the mirror and the cavity field appear to be separable, in 
favour of the quantum correlations between ancillary field 
and mirror. By increasing 6>, the specular situation is ob- 
tained with the entanglement in the mirror-cavity system 
being eventually restored. At = 7r/4 (and = 37r/4) 
a symmetric one-mode biseparable state is found: the 
entanglement in the two subsystems is the same. 

However, the most interesting effects are at the level 
of the mirror state. In terms of equivalent squeezed ther- 
mal states, the interaction with the ancillary mode and 
its subsequent projection onto the vacuum cools down 
the mirror state. This can be put in correspondence with 
the entanglement distributed in the trimodal state. At 
6> = 0, the mirror and the ancilla are not entangled at the 
beam splitter. The situation goes exactly as in the case 
without ancilla and we expect a behavior which repli- 
cates the one corresponding to a no-measurement case. 
However, as soon as ^ ^ 0, measuring the state of the 
ancilla means, effectively, reducing the temperature of 
the mirror, which reaches a minimum at A ujm and 
= 7r/2, when all the entanglement in the trimodal state 
is in the mirror-ancilla subsystem. It is possible to pro- 
vide an analogous study of the behaviour of the effective 
squeezing in the conditioned state. We mention that the 
squeezing function of the conditioned state will initially 
{6 = 0) coincide with the squeezing function of the un- 
conditioned state and then evolves toward a situation of 
large squeezing at small detuning. The behavior becomes 
exactly the same as in Fig. [7] with = 7r/2. 
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IV. CONCLUSIONS 

We have studied the phase-space behavior of an op- 
tomechanical system comprising a movable mirror cou- 
pled to the field of a driven optical cavity, a device which 
is currently of enormous experimental interest. By mod- 
elling the reduced state of each subsystem, in the pres- 
ence of the relevant noise mechanisms, such as a squeezed 
thermal state, we have shed light onto the optomechan- 
ical dynamics originating from their radiation-pressure 
coupling. The optical and mechanical subsystems be- 
have in qualitatively different ways, the cavity field being 
squeezed regardless of the cavity-pump detuning. The 
information gained through this approach allows the de- 
sign of a strategy for the conditional cooling of the mirror 
state via measurements performed on the cavity field. In 
this context, by postselecting the vacuum state of the 
cavity field, we have been able to generate an effective 
measurement-induced non-linearity sufficient to consid- 
erably modify the dynamics of the mirror. Our study 
is a step forward towards the complete understanding of 
the evolution of an open optomechanical system, whose 
dynamical features can be significantly different from the 
unitary dynamics of its noiseless counterpart flS^ . In the 
quest for the coherent use of massive micro-mechanical 
systems operating at the quantum limit, it is mandatory 
to achieve the full control of their coupling to hght, which 
serves as the ideal ancilla for quantum state engineering, 
information processing and state revelation. By suggest- 
ing conditional processes that are able to effectively ma- 
nipulate the state of an inaccessible mechanical system, 
our work contributes to such an important goal. 
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APPENDIX A 

In this Appendix, we give a technical description of 
the general effect that squeezing has on the covariance 
matrix of a single bosonic system prepared in a Gaussian 
state. Let us consider the single-mode squeezing oper- 
ation S'(se*'^) = exp[|e*'^6''"^ — h.c] of squeezing param- 



eter se*'^ (here is the creation operator of a general 
bosonic mode) [23]. This changes a covariance matrix 
V = Diag['Uo, '^i] (which is general enough to represent, 
under proper conditions, a coherent or a thermal state) 
via the corresponding symplectic transformation 

X _ /cosh 5 — cos sinh s — sin(/?sinh5 \ 
^ I — sin cp sinh s cosh s + cos (p sinh s J 

(A-l) 

as {s^(f)vS{s, (f). For cp = or tt, vs (corresponding 
to squeezing and antisqueezing of the bosonic mode, re- 
spectively) reduces to a diagonal form. Therefore, a real 
squeezing factor (either positive or negative) corresponds 
to a diagonal covariance matrix of the squeezed bosonic 
mode. This is the case for the mirror state, whose covari- 
ance matrix M is associated with = tt (cfr. Sec. [TT|) . 

APPENDIX B 

We now give the detailed expressions of the squeez- 
ing parameter and effective mean photon number of the 
mirror state. As is diagonal, it is straightforward to 
prove that e^^^ = uJmE^^/E^ with 

= G^uOmTn/^ + 2Nn{-f + n)5'^[uj^^ - 2ujl^{-fn + A^) 

+ 2-fKS^]} + 2Nf^{-f + f^)[uj^ - 2ujl^{-ff^ + A') + z^S^] 

(B-1) 

with N — 2nM + 1, F = 7^ -h 2/^, 7 = — F, z = 
Y^7^ + A^ and 8 = \/ ti^ + A^. On the other hand, by 
solving (nM + 1/2) e^*^ = Mn using the above ex- 
pression for sm it is easy to obtain the dependence of 
the mean phonon num ber on the detuning. We get 
Hm = -1/2 + {A^/A^)^/ujmE^/E^ with 

= -G^ujmTi^A - 2Ni^S^{j + i^)[uj^ - 2cj^(7/^ + A^) 
+ z^S^] + G^{2Nuj1^kA{j + a^) + 2kS\kz^ - oul^j) 
- NujmA{j + k)[{2j - r)^^ + {2^ + F)A2]}, 
A^ = 2{-G^A + ujmS^)[-2ujtK{j + ,^) + G^uJmT^A 
+ 4.ujl^f^{-f + f^){-ff<i + A^) - 2/^(7 + f^)z^S^]. 

(B-2) 

These equations are operatively very convenient. Their 
exphcit plug-and-play nature allows to quantify the ex- 
pected degree of squeezing and the quantitative thermal 
character of the mirror state per assigned experimental 
configuration. 
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